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' sume that the magnetic field is generated by a causal process, namely a first order phase 

transition in the early universe. As for the topological defects case, we apply the general 
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■r^-j- . relativistic junction conditions to match the perturbation variables before and after the 

, phase transition which generates the magnetic field, in such a way that the total energy C !) 

momentum tensor is conserved across the transition and Einstein's equations are satisfied. 
We further solve the evolution equations for the metric and fiuid perturbations at large 
I scales analytically including neutrinos, and derive the magnetic Sachs Wolfe effect. We find 

^ • that the relevant contribution to the magnetic Sachs Wolfe effect comes from the metric 

^ perturbations at next-to- leading order in the large scale limit. The leading order term 

is in fact strongly suppressed due to the presence of free-streaming neutrinos. We derive 
the neutrino compensation effect dynamically and confirm that the magnetic Sachs Wolfe 
spectrum from a causal magnetic field behaves as £{i -\- 1) Cf oc i"^ as found in the latest 
numerical analyses. 
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1. Introduction 

The origin of the large scale magnetic fields observed in galaxies and clusters is still un- 
known: one of the possible explanations is that they have been generated in the primordial 
universe. A primordial magnetic field of the order of the nanoGauss could leave a de- 
tectable imprint in the cosmic microwave background (CMB) anisotropics. This has been 
analysed in several works: for the scalar mode, see for example |l|, ||, ^, ^, 0]. Here 
we concentrate on the effect that a primordial magnetic field can have on the temperature 
CMB spectrum at large scales: in particular, we focus on the Sachs Wolfe effect. The mo- 
tivation is, that conflicting results are present in the literature regarding the ^-dependence 
of the Sachs Wolfe effect induced at large scales by a primordial magnetic field generated 
by a phase transition: the analytical analysis of ^ found i{£+l) Cf scaling as or more 
negative, and the same result was found in the numerical calculation of Q; on the other 
hand, |, |, |, 0, g found £{£ + 1) Cf scafing as i"^. 
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The aim of this paper is to solve this discrepancy analyticaUy, and we find that the 
relevant contribution to the Sachs Wolfe effect is the one found in |Q, ^, |6|, 0, ^. Naively, 
the magnetic field anisotropic stress would induce a CMB spectrum £{l + 1) oc at 
large scales, basically due to the fact that the metric perturbation $ is proportional to <I> oc 
nB('H/A;)^ at leading order in the large scale expansion k/T-L <^ 1, where IIb is the magnetic 
field anisotropic stress. However, as soon as neutrinos decouple and start to free-stream, 
they develop a non-zero anisotropic stress which adjusts to compensate the one coming 
from the magnetic field, see [^. We demonstrate here that this compensation cancels the 
leading order contribution to the CMB spectrum, and the dominant contribution becomes 
the one from the next-to-leading order in the k/Ti <^ 1 expansion, namely {k/Ti)^, which 
induces then £{£ + l)Cf oc A residual contribution to the Sachs Wolfe remains: 
this is the relic effect of the period of time after the magnetic field generation but before 
the decoupling of neutrinos. However, this contribution is completely unobservable in the 
CMB today. 

In this paper we concentrate on causal magnetic fields, generated at a phase transition 
in the radiation dominated universe before neutrino decoupling (see for example |^0|). We 
make this choice in order to be able to consistently set the initial conditions (in a following 



work we will analyse the infiationary generation case |11|). We assume that the magnetic 
energy momentum tensor is first order in perturbation theory, and that the magnetic field 
has no background contribution: it is a stochastic primordial magnetic field with spectral 
index n > 2 because of its causal generation |jl2| 



{Bi{k)B*{q)) = (27r)^5(k - q)(5,, - kik,)Ak^ , (1.1) 

where A is the amplitude of the spectrum. The magnetic energy momentum tensor is 
therefore automatically gauge invariant. Moreover, we work under the one-fiuid MHD 
approximation, meaning that the conductivity of the universe is high, so that we can 
neglect the electric field and charge separation phenomena, occurring at very small scales. 
Since we are ultimately interested in the CMB spectra at large scales, we also neglect the 
presence of baryons and work under the tight coupling approximation. The universe is 
therefore composed only of a radiation component (photons and massless neutrinos - we 
also neglect neutrino masses), a pressure- less matter component (cold dark matter), and 
the magnetic field. 

In the first part of the paper, we solve Einstein's equations in the long wavelength 
limit, i.e. consistently neglecting all terms proportional to {k/T-Cf'. By doing so, we can 
find analytic solutions for the gauge invariant variables describing the scalar metric pertur- 
bations and those of the total (radiation plus matter) fluid, which hold at scales larger than 
the horizon. Indeed, to calculate the Sachs Wolfe effect, we need the metric perturbations 
at scales larger than the horizon at recombination. 

The solutions are completely determined once the initial conditions have been specified. 
To set the initial conditions, we follow what has been done for the topological defect case 



P^ , 14|. First of all, we assume a sudden phase transition, as for example the electroweak 
phase transition. The magnetic field is then instantaneously generated, if one considers 
perturbations with length scale relevant for the CMB today: k/%-Q <C 1, where T-Lb is the 



Hubble time during the phase transition, corresponding approximatively to its duration. 
Before the phase transition, the universe is the usual perturbed FLRW universe with only 
the standard inflationary adiabatic perturbations; after the phase transition, the magnetic 
field energy momentum tensor contributes at first order to the metric perturbations. In 
order to connect these two stages, we match the geometry and the fluid variables on the 
surface of constant density, so that the induced three metric and the extrinsic curvature 
are continuous [l^ . This implies the conservation of the total energy momentum tensor at 
the magnetic field generation time. 

With this procedure, we find the solutions for the metric and fiuid perturbation vari- 
ables. Each variable is given by the sum of the adiabatic inflationary mode plus the 
magnetic field contribution, which is null before the transition. We have no freedom on 
how to add the magnetic field contribution: this is entirely specified by the matching. 
The solutions are such that the total fluid energy density and momentum are conserved 
at large scales, and the magnetic field does not alter the total curvature perturbation at 
large scales: this is given only by the infiationary contribution. Therefore, the matching 
procedure selects an isocurvature magnetic mode which leaves the curvature unchanged. 

After neutrino decoupling, the solutions change again due to the non-zero neutrino 
anisotropic stress. In order to include analytically the free-streaming neutrinos in our anal- 
ysis, we solve the Bardeen equation combined with the neutrino conservation equations, 
and find an analytical fit for the time evolution of the neutrino anisotropic stress. Then we 
introduce this fit back into Einstein's equation and solve for the metric and fiuid perturba- 
tions including the contribution from the free-streaming neutrinos. With these complete, 
analytical solutions we then evaluate the Sachs Wolfe effect induced by the presence of the 
causal magnetic field. 

We confirm that the free-streaming neutrinos have a fundamental impact on the mag- 
netic Sachs Wolfe |5|, |6|, |^. Neglecting their presence one would conclude that the 
temperature anisotropy due to the magnetic field is proportional to Hb / {kr]i)'^ , where rji 
denotes approximatively the conformal time at recombination. Therefore, the magnetic 
anisotropic stress would completely dominate the large scale Sachs Wolfe effect with re- 
spect to the magnetic energy density, and would induce a CMB spectrum of the type 
£{£ + !) Cf oc £~^. However, as already pointed out in |9|, the neutrino anisotropic stress 
acts to compensate and reduce the magnetic field one. As a result, the Sachs Wolfe con- 
tribution correctly accounting for neutrinos becomes of the form f{r]B,i]v,Vrec)^B/{krii)'^, 
where f{'i]B,'ilu,'nrec) is a function of conformal time at recombination T^reo of the magnetic 
field generation time 77B and of the neutrino decoupling time 77^ , which strongly suppresses 
the magnetic Sachs Wolfe contribution. Therefore, the contribution to the magnetic Sachs 
Wolfe behaving as l/{kr]i)'^ is only due to the evolution of the metric perturbations between 
the magnetic generation time and neutrino decoupling. 

This contribution has been neglected in the literature |^, |5|, |6|, ^ , because it has been 
interpreted as connected to a decaying mode. However, with our analytical approach, we 
can show that the magnetic field solution does not behave as the standard inflationary 
one: one cannot really distinguish a growing and decaying mode, since all modes have 
a comparable amplitude at horizon crossing, and would therefore leave the same imprint 



on the CMB. The reason why this contribution is neghgible only resides in the action of 
free-streaming neutrinos, the anisotropic stress of which counteracts the magnetic field one 
and effectively cancels the large scale temperature anisotropy. The new contribution to the 
Sachs Wolfe effect can only act as a source of anisotropy after the magnetic field generation 
time but before neutrino decoupling, and is negligible. Therefore, we finally confirm the 
result present in the literature [Q, |5|, ^, |^] : the dominant contribution to the magnetic Sachs 
Wolfe effect becomes the one coming from the metric perturbations at the following order 
in the long wavelength expansion, i.e. at order {krji)^ . This is what has been calculated in 
all the numerical analyses of the CMB temperature anisotropy. In the second part of the 
paper, we therefore proceed to evaluate analytically this contribution. 

In order to solve Einstein's equation at next-to- leading order, i.e. including terms of 
the order {k/T-Cf , the easiest way is to solve for the curvature perturbation equation. This 
allows us to find the solution for the Bardeen potentials at order (kiji)^, and consequently 
all the other fluid variables. We can then evaluate the Sachs Wolfe contribution at this 
order. 

The structure of the paper is the following: in section ^ we derive Einstein's and 
conservation equations in the presence of a magnetic fleld. In section ^, we explain the 
matching procedure, used then in section ^ to calculate the metric and fluid variables 
at leading order in the kr]i <C 1 expansion, both before and after neutrino decoupling. 
In section ^ we compute the Sachs Wolfe effect from the leading order solution, and we 
conclude that the relevant contribution comes from the next-to-leading order. Therefore, in 
section ^ we calculate the next-to- leading order solutions, and we conclude in section 0. We 
consider scalar metric perturbations on a spatially flat Friedmann background, and we work 
with gauge invariant variables using the notations of ||l^ (with respect to [jl^, we have $ 
with the opposite sign while ^ is the same). We normalise the scale factor to one at equality, 
so that in a matter plus radiation universe a = y'^ +2y in terms of the dimensionless variable 
y = v/Vii where ry denotes conformal time, rji = ?7cq/(\/2 — 1) ~ rjj-cc and r^cq represents 
conformal time at equality. Scales larger than the horizon at recombination satisfy xi = 
krji ^ 1. A dot denotes derivative with respect to conformal time, while a prime denotes 
derivative with respect to y. Greek indexes go from to 3, while latin ones from 1 to 3. 
The neutrino background energy density fraction is denoted = Pu/Prad — 0.4. 

2. Einstein's and conservation equations 

In this section we derive Einstein's and conservation equations with a non-zero magnetic 
field. We assume that the electromagnetic energy momentum tensor is first order in per- 
turbation theory, meaning that the fields themselves are half order. The electromagnetic 
field tensor and the energy momentum tensor can therefore be defined with respect to the 
unperturbed velocity of the fiuid energy frame of the FLRW background, m° 
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where e^^^ = S^^^fj^^] is the totally antisymmetric rank 3 tensor. The components of the 
energy momentum tensor are 

p^^ = ^{E^ + B^) , (2.4) 

Pem = UE^ + B^), (2.5) 
b 

q^^ = e'^p^E^B^ , (2.6) 



a 



+ 5 - E'^Ep - B^Bp . (2.7) 
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From Maxwell's equations F'^^^/s = j" one sees that the current density j° is also half 
order as the electromagnetic field. The covariant Ohm's law can therefore be written with 
respect to the background velocity 

r + n°n/3/ = aF"% . (2.8) 

We work under the MHD approximation, i.e. the conductivity a is infinite. Prom the above 
equation one sees that, in order to keep the current finite, the electric field must vanish 
1 17]. Therefore in the following we set the electric field to zero. The MHD approximation 
is valid on sufficiently large scales, where charge separation effects are not important. 

The total energy momentum tensor includes the fluid, labelled by f and representing 
radiation and matter, and the magnetic field. Using the notation p^ = B'^/2, = 
TTg • = {B'^/3)g^j —B^Bj, the components of the total energy momentum tensor in real space 



3j 
are 



r% = -PF - SpF - PB , (2.9) 
T'o = -a{pF+PF)5u\ (2.10) 
T'j = {pF + SpF + Pb)4 + PF ttf'j + ttb'j ■ (2.11) 



We only consider scalar perturbations. Following |T^, we expand the scalar part of the 
metric and fluid perturbations by scalar harmonic functions Y = e"'*' -^ (for the deflnition 
of y* and Y^j see |jl5[). The perturbed scalar energy momentum tensor in wave- number 
space becomes then 

6T% = {-6pF - pb)Y , (2.12) 
6To = -{pf+pf)vY\ (2.13) 

6Tj = {5pF + PB)5i + {pF TTF + TrB)Y'j , (2.14) 

where v is the scalar part of the velocity perturbation (5n* and (see [|l|, [l8| ) 

p^ = ^{Bi*B'), (2.15) 
= lk%{B, * B^) - ^(i?„ * B'n , (2.16) 



The scalar perturbed Einstein equations 5G^iy = SnGST^i, in terms of gauge invariant 
variables, in the presence of the primordial magnetic field and in a spatially flat universe 
i m are 

A:2$ = AnGa'^pFiD + — ) , (2.18) 

PF 

-A;2($ + ^) = SirGa'^pFiTTF + — ) , (2.19) 

PF 

k{n^-^) = AttGo^ {pf+Pf)V , (2.20) 

where $ and ^ are the Bardeen potentials, D is the gauge invariant variable corresponding 
to the density perturbation in the velocity-orthogonal slicing, and V is the gauge invariant 
variable corresponding to the velocity perturbation in the Newtonian longitudinal gauge 



The conservation equations Tp ^ + Tg ^ a ~ ^ involve the momentum exchange Qi 
between the fluid and the magnetic field, represented by the Lorentz force (the energy 
exchange is set to zero since it involves the electric field). Projecting the conservation 
equation on the rest space of the comoving observer by means of the projector /i"^ = 
9^/3 + gives 

ha^TF} = -hjT^l = e^^sJ^B^ = Qa , (2.21) 

where Qa is a purely spatial vector (Qq = 0) denoting the Lorentz force. We Fourier 
transform it and extract its scalar part by —ik'^Qi = £bY, which from Eq. (2.21) must 
satisfy the relation 

^ = ^B + -%, (2.22) 
£b = ^{Bi* B') - k^kiiBj * B') . (2.23) 
The fluid conservation equations become 

O-TJ 

b - swnD = -k{i + w)v - 2wn'KF + t^Ib , (2.24) 

kpF 

V + nV = k-f + -^D-l--^k7TF + —^^, (2.25) 

1 + W 31 + W PF{1 + W) 

where w = Pf/pf, = Pf/ 'Pf and the fluid has no internal entropy perturbation Ep = 
(we remind that we have set the energy exchange between the fluid and the magnetic field 
to zero). We introduce the following notations 

^B = ^ Ub = ^ Lb = ^, (2.26) 

Prad Prad Prad 



with pf = Prad + Pmati SO that Eq. ( p. 22 ) becomes Qb = 2nB + 3Lb//c. 
3. Matching conditions 



We assume that the magnetic field is generated by a causal process acting 'fast', i.e. within 
one Hubble time, as for example a sudden phase transition in the early universe. We denote 



the magnetic field generation time by r]B- Before this time, only adiabatic perturbations of 
inflationary origin are present in the universe. Afterwards, there is an extra contribution 
to the metric perturbations due to the magnetic field, which in turns affects the fiuid 
perturbations. In order to establish the initial conditions for this system, we follow what 
has been done in |13| for the analogous case of the topological defects. They perform 
a matching of the perturbations on a constant energy density surface, which allows to 
make the link between the pre- and post-magnetic field phases in the correct way: the 
physical configuration satisfies the conservation of the total energy momentum tensor (and 
consequently Einstein's equations) at all times. The two independent matching conditions 
are given by Eqs. (36) of [|l3|, and read 



where 



and 



[SnGal ipFDg + pB)]^ = 



F± = lim [F(r/B + e) - F(r/B - e) 



Dg = D - 3{1 + w)—V + 3(1 + w)^ 
K 



In terms of the variable D and using Eq. (|2.18D we obtain then the conditions 



Z)(77B + e) + 



D{iiB - e) 



1 + ae 
V{m + e) = V{iij, - e) . 



(3.1) 
(3.2) 

(3.3) 
(3.4) 

(3.5) 
(3.6) 



We now proceed to solve the system of Einstein's plus conservation equations at large 
scales. We need to split the universe evolution in three stages: before magnetic field gen- 
eration, after magnetic field generation and after neutrino decoupling. For the transition 
from the first to the second stage we impose the matching conditions derived here, while 
neutrino decoupling does not introduce any discontinuity in the metric and fiuid perturba- 
tion variables. 

4. Leading order solutions at large scales xi ^ 1 

In order to solve the system of Einstein's plus conservation equations, we choose to derive 
a second order differential equation for the variable D. We combine the total density and 
momentum conservation equations (2.24) and (2.25) to get 



£) + (1 + 3cf - Qw)nD + 2,%^ 

w / k 



1 



--Aw + -w^ + 3c: + -f — 



271 



-2w + 3ci + 3w'^ + 



3 \n 



k 



n 



D 



TTp — 2'Hw'k-p + 



l + a 



3c2 + 



a 1 / k 
l + a '^3\n 



(4.1) 



In order to determine the Sachs Wolfe effect, we only need to solve for scales which are 
over the horizon at recombination, xi = ktji <^ 1. Therefore, in a first instance, we drop 
the terms proportional to {k/T-Cf'. 

4.1 Before neutrino decoupling 

Before their decoupling at a temperature of Tjy ~ 1 MeV, neutrinos do not free-stream and 
have therefore zero anisotropic stress. They are characterized only by their background 
energy density pu, their density perturbation Dy and their velocity perturbation Vy. These 
are simply included in the total matter perturbations D and V . Consequently, for T > 
the fluid anisotropic stress vrp in Eq. (|4.1| ) vanishes. 

In order to solve Eq. (^T|) in a matter plus radiation universe with = Prad + Pmat , we 



follow I IE]. Before magnetic field generation the source in Eq. ( [4.1[) drops, and we have the 
usual homogeneous solution, that we express here as a function of the variable y = r]/r]i 



D-{y) 
usiy) 



aiURiy) + a2Us{y) 
1 

2y + 3y2 + y3 ' 



y3(8 + 3y)(10 + 10y + 3y2 



(4.2) 
(4.3) 

(4.4) 



9(1 + y) 

where ai and 02 are arbitrary constants and we set the decaying mode to zero: 02 = 0. 
We remind here that in terms of y the scale factor is a = + 2y and therefore the 
Hubble parameter is H'^ = 4(1 + a)/{rjla^) = 4(1 + y)'^ /[rjfiy'^ + 2j/)^], the equation of 
state parameter of the total fluid is w{y) = 1/[3(1 + y)^], and the total sound speed is 
c2(y) = 4/[3(4 + 6y + 3y2)]. 

In order to solve after the generation of the magnetic field we use the Wronskian 
method, where the source is given by the right hand side of Eq. (4J) with ttf = in terms 
of y 



1 + y)^ (1 + yf{8 + 6y + 3y^)nB + 2(8 + 26y + 37y^ + 24y^ + 6y^)UB 



The solution becomes then 



y2(2 + y)2(4 + 6y + 3y2) 



D+iy) 
W{y) 

Bi{y) 

B2{y) 



[61 + Bi{y)]uniy) + [62 + B2{y)]usiy) , 



where 



dus du^ 

dy dy 

, usSb 
dy 



dy 



W 
urSb 
W 



(4.5) 

(4.6) 
(4.7) 

(4.8) 

(4.9) 



where again bi and 62 are arbitrary constants. The metric perturbations ^ and ^, and the 
velocity perturbation V can now be calculated using Einstein's equations ( 2.18| ) to ( 2.2C| ) 
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both before the magnetic field generation, setting p-Q = ttb = and using given in 
Eq. ( [4.2[ ), and after the magnetic field generation, using instead Eq. (^). Consequently, 
<I>, ^ and V after the magnetic field generation are functions of bi and 62- We then use 
the matching conditions Eq. ( |3.5D and ( |3.6| ) to determine the free parameters 61 and 62. 
This completely specifies the magnetic contribution to the metric and fluid perturbations 
without ambiguity. The only free parameter remains ai, the amplitude of the scalar met- 
ric perturbation from inflation, which is then determined by the COBE normalisation. 
Denoting j/b the magnetic field generation time, the matching gives 



3 16 
ai--nB and ^2 = ^^f^B + 4(1 + yB)nB , (4.10) 



so that 



WTWTW " (TT^ " y(2 + y)(i + y)"'^ • 

The first part of the above expression, proportional to ai, is the inflationary mode of 
Eq. ( [4. 21) , while the second part proportional to JIb and IIb is the magnetic field con- 
tribution. For later convenience, the magnetic field contribution is denoted D^, so that 
D+{y) = D-{y)+D^{y). Note that D^{y) is a decaying mode, but its amplitude is of order 
Ob + 2nB as long as y < 1, i.e. anytime before recombination (in the conclusion section |^, 
we discuss the role of the decaying mode of the perturbations generated by the magnetic 
field). The metric perturbations $ and and the velocity V follow from (|2.1S - 2l20| ), and 
are given in appendix 

Note that these solutions hold at large scales k <C This is because in Eq. ( [4.1] ) 

we have dropped the terms proportional to {k/H)'^ with respect to those constant in k. A 
consequence of this is, in particular, that the momentum conservation equation ( p.25| ) is 
only satisfied at leading order in xi = krji: i.e. neglecting the Lorentz force which is order 
0{x^) and the term proportional to kc^D which is order 0{x\), and keeping only the term 
k"^ which is order 0{\/xi). 

4.2 After neutrino decoupling 

After decoupling, neutrinos start to free-stream, and acquire a non-zero anisotropic stress 
TTj/. Consequently, for r] > rj,^ Eq. ( [4.11 ) gets a new source term, generated by vrp = RuT^u 



(where Ru = Pu/Pt&a)- Since, however, the neutrinos do not add as a new component in 
the universe (as did the magnetic field), in this case we do not need to perform a matching 
to guarantee the conservation of the total energy momentum tensor. The new anisotropic 
stress component does not introduce a discontinuity, since it builds up at sub-horizon scales 
continuously in time, so that iru{yu) = 0. We have then 

D~{y) for y <yB , 

D{y) = { D+iy) = D-{y) + D^{y) for ye < y < y. , (4.12) 



where D (y) is the inflationary solution Eq. (|4.2[ ) with 02 = 0, D^{y) is given in Eq. ( 4.11J ) 



and D^{y) is generated by the neutrino anisotropic stress, such that D^^{y^) = 0. D^{y) 



is calculated again using the Wronskian method, where now the source is given by the fluid 
anisotropic stress in Eq. 



Su{y) = -'^{niiif \-2w + 2,cl + 3u;2l i?^7r^(y) + 2nr]iwRu^ ■ (4.13) 
L J dy 

We again consistently neglect the terms proportional to {k/'H)'^. The solution becomes 
then 

D'^^y) = Ci{y)uRiy) + C2{y)usiy) , where (4.14) 

As explained above, since the source Su{y) is continuous, the energy-momentum conser- 
vation is automatically satisfied at the neutrino decoupling time y^, and no additional 
matching conditions are needed. Hence once the neutrino anisotropic stress T^u{y) is deter- 
mined, the solution D^^ is completely fixed. 

The evolution of iiy in the presence of an external constant anisotropic stress has been 
studied in [^] (see also [^). Even though the neutrinos and the magnetic field do not 
interact directly but only through gravity, the neutrino anisotropic stress quickly adjusts 
to the external one and compensates it. In order to determine the time evolution of Hy, 
we follow the method given in j^, i.e. we combine the neutrino conservation equations and 
Einstein's equations to derive a fourth order differential equation for iTu{y)- This derivation 
is presented in appendix The final time dependence of 'Ku{y) is rather complicated, but 
it can be approximated well by (c/. Fig. Q) 

SHb (yl \ 40 ai 

With this approximation, we can solve analytically the integrals in Eq. ( |4.15| ). We find 
then 

nas/,,x 80i2y + y,)iy - 4(y - y,f 

^ - 9(15 + 4W + ,)(2 + ,)^^ + yM + y){2 + y) " ^''''^ 

The part proportional to ai is the standard neutrino contribution to the inflationary per- 
turbation: we denote it Df^^. The part proportional to IIb is the contribution due to the 
interplay among the magnetic and neutrino anisotropic stresses: we denote it D^^^, so 
that D^{y) = D^^{y) + -Dmag(y)- The metric perturbations and and the velocity 
V^^ are again given in appendix It is clear that the time dependence of Eq. ( 4.17] ), and 



consequently of the other metric and fluid perturbations given in the appendix, follows 
directly from our choice for the approximated form of TTy{y) given in Eq. ( [4.16 ). Since this 



fit interpolates well the behaviour of the real neutrino anisotropic stress, these expressions 
will give the right order of magnitude for the evaluation of the Sachs Wolfe effect. However, 
their detailed time dependence is to be considered only indicative. 

Summarising, we have that the final metric and fluid perturbations valid after neu- 
trino decoupling can be decomposed into three contributions: the standard inflationary 
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Figure 1: Blue, solid line: the neutrino anisotropic stress as a function of y derived in appendix 
Eq. ( B.12 ), and red, dashed line: the fit given in Eq. ( 4.16 ). We have chosen comparable values 
for Ub 
at y : 



oi, and normalised by IIb. The two curves flatten at 7ry(j/) — —SYIb/Ru- The growth 
1 represents the inflationary contribution to the neutrino anisotropic stress, i.e. the second 



term in the right hand side of Eq. (4.16). The time at which this contribution starts to dominate 
depends on the relative amplitude of IIb and ai. 

one, the magnetic one arising after magnetic field generation, and the one due to the 
free-streaming neutrinos which further contribute after their decoupling to the metric per- 
turbations through their anisotropic stress 



D^'^iy > y. 



D-{y)+D^{y)+D'^{y) , 
<I>-(y) + <I>i^(y) + <i>-(y) , 
vI/-(y) + ^B(y) + vI/-(y) , 
V-{y) + V''{y) + V^{y) . 



(4.18) 
(4.19) 
(4.20) 
(4.21) 
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o 
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o 
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Using the above solutions we can calculate the curvature perturbation By imposing 
the continuity of the induced three metric and the extrinsic curvature on the surface of 
constant density, the matching conditions guarantee in particular that C is continuous at 
the magnetic field generation time ijb- Consequently, the magnetic contribution to C is 
by definition such that it vanishes at ijb, and only the constant inflationary part remains. 
Moreover, if we calculate C for rj > rj-Q using the solutions for <I> and ^ derived above, we 
further find that the magnetic contribution to it vanishes at all times 

This shows that the curvature perturbation is conserved at large scales xi <C 1: therefore, 
if the magnetic contribution is zero at rye, it must remain so for rj > rj-Q aX, leading order in 
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xi, i.e. C oc 0{l/x\). The neutrinos also do not contribute to the curvature perturbation, 
since we have imposed adiabatic initial conditions for the inflationary solution 



C,{y > Vu 



3(1 + ^1-) 



fin 



fin 



10 ai 



n 



(4.23) 



We will see in the following that the curvature is no longer conserved for intermediate 
scales, i.e. at the following order in the xi <C 1 expansion, C oc 0(2;?) (c/. section 1). 



5. Sachs Wolfe effect from the leading order solution 

In the approximation of instantaneous recombination, valid for wavelengths larger than the 
Hubble scale at recombination xi ^ 1, the Sachs Wolfe contribution to the temperature 
anisotropy is (see e.g. p9[|) 



AT Dy^ 

— {k,7]o) ~ ^ 



(A;, 



Tec J ■ 



(5.1) 



In order to evaluate Dg^, knowing the Bardeen potentials $ and ^I* which have been derived 
above, we use the conservation equations for the photon fluid 

4, 



97 



= - $) + -Dg^ + 



-Lb 



(5.2) 
(5.3) 



These follow directly from (2.24) and ( |2.25 ). We remind that we work under the tight 
coupling approximation, so that we neglect the Thomson scattering term, we equal the 
baryons and photon velocities ~ Vb, and neglect the baryons energy density so that the 
parameter i?b = 3/9b/4/?-y = 0. Note that in tight coupling the momentum exchange due to 
the Lorentz force must be included in the photon momentum conservation equation, since 
it acts on baryons and ~ Via. From the above equations it follows that the Sachs Wolfe 
contribution is 

AT ^ Dg -yiVr 

T ~ 4 



+ *(yr, 



$(2/rec) 



4(1 -i?.) 



Lb 



(5.4) 



we therefore only need to evaluate Vy. In order to do this, we differentiate Eq. (|5.3| ) and 
substitute with Eq. ( ^ ) (note that Lb = ^(Ob — 2nB)/3 is constant) to obtain a second 
order differential equation for V^, which reads in terms of the variable y 



dy^ + 3 



d 



$) . 



(5.5) 



Following what found in the previous sections and presented in the end of section |4.2| , it 
appears that the source term can be divided into the standard inflationary contribution. 



the magnetic contribution, and the free-streaming neutrino contribution (c/. Eqs. ( A.l 
A. 2 ) and (A^, [A.5| ) of appendix^. Consequently, can also be divided as 



o 
o 

o 
I— ^ 
o 

o 
o 
o 



K = v: 



(5.6) 
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The standard inflationary contribution can be determined by solving directly the differ- 
ential equation ( |5.5| ), where the right hand side is given by ^'"(y) and ^~{y) {cf. Eqs. ( [A.l[ ) 
and (A. 2) of appendix The magnetic and free-streaming neutrino contributions can be 
determined instead using the Wronskian method, and after one integration by parts of the 
source integrals they can be written as 



xi cos 



+ xi sin 



xiy\ 

xiy 
V3 



dy'i^^ - cos 



Vb 



xiy 
V3 



dy\^^ - sin 



Vb 



xi cos 



/ Xiy 
\V3 



+ xi sin 



dy'i"^' 



$3 



cos 



$3 



I sm 



xiy 
V3 



xiy 

xiy' 
V3 



(5.7) 



Note that the boundary of the integrals in Eqs. (|5.7D and (^) are a consequence of the 
fact that Yy is continuous both at magnetic field generation time y^, as implied by the 
matching conditions, and at the neutrino free-streaming time yu. 



As an example we first evaluate the standard inflationary contribution to the Sachs Wolfe 
effect, given by V~. We impose adiabatic initial conditions, such that at large scales xi <^ 1 
we have Vy{y) ~ Vu{y) — ^mat(y) — ^"(y); where V~ denotes the total velocity perturba- 
tion from inflationary perturbations given in Eq. ( [A.3| ). The initial conditions completely 
specify the arbitrary constants ci and C2 coming from the integration of Eq. (5.5). To 
the standard inflationary contribution we should further add the free-streaming neutrino 



contribution, represented by the part of Eq. (5^) sourced by ^f^^ — ^f^f. This is usually 



neglected at large scales, since it has an impact only of a few percent on the standard 
inflationary Sachs Wolfe: we evaluate it in appendix From Eq. (5.4), the inflationary 
Sachs Wolfe effect is then given by 



AT 



inf 



Using Eq. (|A.3[) we find 
AT 

Comparing with Eq. 



inf 



2(40 + 40j/,ec + 18yl^ + Syj^) ai ^ 2.5ai 
3(2 + yrec)3 xj"^ xf ■ 

one sees that this result corresponds to the usual $~/3 



(5.9) 



(5.10) 
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We now calculate the Sachs Wolfe effect due to a non-zero primordial magnetic field. 
Combining Eqs. (|5.4|) and (|5.6|) this reads, in terms of the variable y 



AT 
~T~ 



1 

xi 



dV} 



+ 



dV^ 



dy dy 



3 Lb 
4(1 -i?^) k 



(5.11) 
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where for V^'^ we only take into account the magnetic (V^) and not the inflationary part 
(VJ^£, calculated in appendix Q). Differentiating Eqs. ( |5.7] ) and ( |5.8D and further writing 
the Lorentz force as L^/k = f^e/S + 2113/3 we finally find 



AT 



+ 



V3 



cos 



sm 



xiy 
^/3 
xiy 
V3 



as 

mae 



1 



4(1 - 



■(J7B-2nB) 



(5.12) 



dy'(*^ 
dy'(*^ 



sin 



cos 



xiy' 
V3 



+ 



+ 



y 



dy'i^ 



mag 



mag 



mag 



mag 



) sin 
) cos 



y-^ 



xiy 

xiy' 
V3 



The leading term at large scales xi ^ 1 comes from the first two terms of the above 



equation, which can be derived from Eqs. (A.l) and (A. 2). We immediately see that it is 
of order ©(l/xf), while both the Lorentz force term and the integrals are of order 0{xi). 
Let us then concentrate on this apparently leading term. For the purely magnetic part, we 
find 



12 i^y' +2y- 4(1 + y)yB] Hb 
y3(2 + y)3 xj 



(5.13) 



Supposing that the magnetic field is generated at the EW phase transition at T ~ 100 
GeV, one has ys = "Hb/iIi — 10~^^, while yrec — 1- Therefore, the above expression does 
not depend strongly on the magnetic field generation time. For the neutrino 'magnetic' 
part, we find instead 



mag 



12 [(3y2 + 2y - (6 + 5y)y,){y - y,)] 



mag 



(5.14) 



y4(2 + y)3 xj ■ 

Here we can see at play the compensating effect of neutrinos already proposed in [§, |5|, |9| 
(see Fig. ||) . The contribution in Eq. ( 5.13| ) that is not suppressed by yB is exactly cancelled 
by the effect of the neutrino anisotropic stress. For the total Sachs Wolfe effect at leading 
order in the xi <C 1 expansion 0{x^'^), we find then 



AT 



12 



4yrcc(yB - 2y^) + 2/rec(4yij + 5yl - 8y^) 



Hb 



yrec(2 + 2/rec)' 



XT 



(5.15) 



This contribution is proportional to the time of generation of the magnetic field ys — 10~^^ 
and to the time of neutrino decoupling, for which we have y^ ~ 10~^: it is therefore 
strongly suppressed. It actually only corresponds to the imprint of the magnetic field 
anisotropic stress from its time of generation yB to the decoupling time of the neutrinos 
y^. The subsequent magnetic contribution to the Sachs Wolfe effect, arising from y^, up to 
recombination time, is cancelled by the free-streaming neutrinos. 

The dependence on y and y,^ of Eq. ( |5.15 ) follows from the particular form of the 
potentials (A.4) and ( [A.5D , which in turns depends on the fit function TT,/{y) given in 
Eq. (4.16). Using this fit and setting y^-ec = 1 gives the value 



AT 



-7-10 



-6 Hb 



(5.16) 
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This becomes of the order IIb only at scales k < 10^^ Mpc^^: larger than the present 
horizon! A different fit could provide a different numerical value but would not alter the 
overall suppression effect due to the free-streaming neutrinos. For example, if we had 
assumed that iru{y) adjusts itself linearly to the compensating plateau —SUb/Ru, instead 
of quadratically as chosen in Eq. (4.16), the numerical value of the temperature anisotropy 
( 5.16 ) would have changed by about 14%. Even though the fit of Eq. ( 4.16| ) is better than 
a linear adjustment, conservatively the precision of our analytical analysis this should be 
taken to be of about 10%. 

Note that, if we neglected neutrinos, the magnetic Sachs Wolfe effect would be given 
by Eq. ( ^.13 ) taken at recombination time, y = y^ec- This gives approximately 



AT 



20 Ub 



(5.17) 



B ^ J-i 

We will see in section ^ that this is of the same order of magnitude of the next-to-leading 
order contribution to the magnetic Sachs Wolfe. It would then have an observable effect 
on the CMB for a magnetic field of the order of the nanoGauss. It is only because of the 
equilibration among the magnetic and neutrino anisotropic stresses that this contribution 
is washed out, and not because it comes from a decaying mode (c/. discussion in the 
conclusion). If we had neglected neutrinos, as is usually done for the infiationary mode at 
large scales, we would still have found the correct amplitude for the magnetic Sachs Wolfe, 
but not the correct spectrum, as we now show. 

5.1 CMB spectrum from the Sachs Wolfe effect at leading order 

With the above result we can give an estimate of the Sachs Wolfe contribution to the CMB 
spectrum. For the scope of this paper, and in order to give an interpretation of our result, 
we are merely interested in its general ^-dependence. We do not aim at determining the 
amplitude of the magnetic contribution given in Eq. ( ^.16 ) to the CMB spectrum at large 
scale in any details, since this is anyway completely unobservable. 
The CMB spectrum is given by (we use the notations of ||20|) 

1 



dkk Qe{rio,k)Q}{r]o,k) , 



^(2^ + 1)2 

where at sufficiently large scales (c/. Eq. (|5.lD and the appendix of l2l[|) one has 



2£ + l 



+ ^{k,r],^ 



<^{k,rir, 



(5.18) 



(5.19) 



where the integrated Sachs Wolfe effect is neglected, and ^(r/rec) denotes the visibility 
function. Since our treatment resides on the hypothesis of a causally created magnetic 
field, we assume that the magnetic and infiationary perturbations are not correlated. The 
magnetic contribution to the temperature anisotropy is then given by Eq. ( |5.15 ): since we 
are interested here only in the ^-dependence of the CMB spectrum, for brevity we resume 
the amplitude of Eq. ( 5.15|) in a generic function /(?/rec, yi/, ye)- We have then 



e,^(7?o,A:) 
2£ + l 



fiyr 



,yi',yB) — 



(5.20) 
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The magnetic CMB spectrum becomes then 



where |nB(A;)p denotes the spectral amphtude of the magnetic field anisotropic stress 

(nB(k)n^(q)) = (2^)35(k - q)|nB(A:)|2 . (5.22) 

This has been calculated in ||l], ^, |l^, and it shares the same /c-dependence as the spectral 
amplitude of the magnetic field energy density J^b- For a causal magnetic field with n > 2 
it is simply constant in k up to the damping scale ko, which we assume time-independent. 
We denote its amplitude by fl, so that 

|nB(fc)|' = n^^-^ . (5.23) 



Substituting in Eq. ( ^.21 ) and setting x = krjQ we find then 

iii + l)Cf - f\yr.c,yu,yB)9HVrec)n^-^^^-^^ dx (5.24) 

^ Pr«H vi {mkor Jo X 



P?ad ^1 imkof Jo 

/ {yicc,yu,yB)g (??rec)n 



where for the last equality we have taken the limit rj^kj:, ^ 1 {cf. [^). Therefore, we find 
a scaling for the CMB spectrum as \/i for ^ > 2, as found in |l|, Note that, although 
this scaling might seem unusual, this contribution is of the same order of the standard 
inflationary one in the xi <C 1 expansion, i.e. ai/xf {cf. Eq. ( ^.lOP ). Contrary to the 



magnetic field case for which the spectrum is flat {cf. Eq. ( 5.23 )), the inflationary generated 
perturbations have an almost Harrison- Zeldovich spectrum corresponding to |ai(fc)p oc k, 
which leads to the usual flat CMB spectrum at large scales once inserted into Eq. (|5.21 ) 
(see Eqs. (3.112) and (3.113) of pl). 



The ^-dependence of the Sachs Wolfe effect coming from Eq. ( 5.15 ) does not corre- 
spond to the result presented in Q |5|, ^, ^: in these works only the next-to-leading order 
contribution, constant in xi, has been considered, leading to a CMB spectrum which be- 
haves as £{£ + i)Cf oc i"^. This is due to the different way of treating the initial conditions. 
Indeed, @, H; only need to derive the initial conditions for the metric and fluid vari- 
ables after neutrino decoupling, which are then inserted into the Boltzmann code. The 
contribution to the temperature anisotropy of the magnetic field anisotropic stress from 
its time of generation y^ to the decoupling time of the neutrinos y^ is therefore absent. 
Moreover, the time evolution of the metric and fluid variables at order ©(l/xf) which we 
have obtained above has been identified with a decaying mode, and therefore neglected 
since the beginning {cf. discussion in section 0). When one derives the initial conditions in 
the syncronous gauge, as done in [Q, ^, 0], it results that the only way to avoid the long 



wavelength mode at and select the constant mode at order 0{xi) is to solve the 

Einstein equations setting the neutrino anisotropic stress to be constant in time, and ex- 
actly equal and opposite to the magnetic field one. This is why in ||, ^, |^ the anisotropic 
stress of the magnetic field is compensated by the one of the neutrinos since the beginning, 
directly in the initial conditions. In our derivation, the anisotropic stress of the neutrinos 
is instead time dependent: we derive its behaviour in time in appendix 1^ and then insert 



it in Einstein's equations by means of the analytical fit given in Eq. (4.16). Therefore, here 
we obtain the compensation effect dynamically, without having to insert it directly in the 
initial conditions. 

Within our treatment, we can somehow mimic the result of [@, |5[ |^, 0] if we set tTi, = 
— SIIbZ-R;/ and y^, = ye^- In this case, the source term ttf = Ri/TTu in Eq. (|4.1| ) exactly 



cancels the magnetic source term He. We can then repeat what done in Sec. [4.1| , where 



the source Sb of Eq. (4.5) is now only given by the l^B-part. The matching conditions are 



different and the solution for D becomes simply 



The source term in Eq. ( p. 12 ) then vanishes at leading order ©(l/xf) and therefore the 



Sachs Wolfe contribution is exactly zero at this order, as found in [pi| ]. 

Our analytical solution ( 5.15| ) shows that the period between the generation of the 



magnetic field and the time of neutrino decoupling would leave an imprint on the Sachs 
Wolfe that is diverging at large scales. However, since this imprint is suppressed by ye and 
yi, as a consequence of the neutrino compensating effect, it leaves no observable impact. 
Therefore, we now proceed to compute analytically the next order contribution 0{x^) in 
the Sachs Wolfe. 



6. Next-to-leading order solutions at large scales xi ^ 1 



The temperature anisotropy at constant order in xi is given by Eq. ( ^.12 ), where we now 



need to take into account the next-to-leading order 0{x^) in — <I>^ and ^'mag ~ *^mag- 
The Lorentz force term, which could be neglected in the above result, must be taken into 
account now, since it is constant in xi. The integrals are one order higher with respect to 
their integrands, so to get the contribution constant in xi it is enough to keep the order 
0{x^'^) in the sources - and ^'^^g - ^mag- 

The easiest way to compute the order 0{x^) in $ and ^' is to use the curvature 
perturbation Indeed, starting from the solution given in Eq. ( [4. IT] ) for D^"^, we can 
calculate the order 0{xj^) in Once this is known, we can use definition (4.22) to integrate 
for $ at the same order. 



""^Note that the condition yi, — j/b is necessary if the neutrino anisotropic stress adjusts itself instan- 
taneously to the magnetic field one. A sudden jump in the neutrinos anisotropic stress induces in fact a 
discontinuity in the total energy momentum tensor, and we can only arrange this by means of including it 
in the matching at i/b- 



Deriving Eq. (|^), using definition ( |2.20D and the momentum conservation equation 



(2.25), we find the evolution equation for ( 



dc _ nm 

dy 1 + w 



2nB 



3(1 + a) 



-WTT-p 



(6.1) 



Since the leading order in D is 0{x^, the source term in the above equation is also of the 
same order: therefore, we can solve the equation to find C, at order constant in xi. 

In solving Eq. ( |6.1D we concentrate here only on the magnetic part and neglect the 
inflationary one, which is not relevant for the scope of the paper. This means that, for 
Vb < y < Uu, the source is given by the variable {cf. Eq. ( 4. 111) ). Furthermore, we 
have vrp = and we solve Eq. ( |6.lD imposing the continuity of ( at ye: for the magnetic 
part this means C^Cj/b) = 0, which is equivalent to the matching conditions. In the limit 
ye < y ^ 1, the solution reads 



IIb 
2 



1 _ ^ 



21og — 



(6.2) 



This result is consistent with what given in Eq. (84) of Q: the curvature is sourced only by 
the magnetic field anisotropic stress. The reason for this is that in Eq. (|6.1| ), initially, the 
fluid energy density is compensated by the magnetic field energy density VL-q because 
of the matching condition (^^) and because ~ 1/3 and ttf = 0. Therefore, the only 



possible source for Q is proportional to IIb- With the evolution, the compensation of D is 
no longer perfect, and becomes proportional to IIb, Eq. (4.11). Therefore, the curvature 
only depends on IIb- 

After neutrino decoupling, we use our approximation ( 4.161 ) for and insert it in 
Eq. (6.1). The rest of the source is now logically given by + -Dmag- before, we 
impose the continuity of Q at using the solution determined previously, C,{yu) = C^{yv)- 
We find, for ?/b < < 2/ < 1, 



C(y^ <y<l) 



0{x\) 



as 

mag 



IIb 
2 



-1 + 2^ 



21og^ 

y \yB 



(6.3) 



which reduces to Eq. (^]^) for y = y^. The curvature stops growing after neutrino decou- 
pling, as already pointed out in since the source of Eq. ( |6.1D is fully compensated: the 
neutrino anisotropic stress compensates the magnetic field one, and -D^ag compensates the 
part proportional to IIb hi which acted as a source for y < y^ {cf. Eqs. (|4.1l|) and 



( 4.17 )). In the limit y ^ 1 > y,y we find instead 



C(y>l>y.)lo(,0) = C'' + C 



B 



mag 



log 



ys 



(6.4) 



This shows that after equality, while the neutrinos are still a relativistic component and 
therefore continue to compensate the magnetic field anisotropic stress, the fluid is no longer 
relativistic, and no compensation is possible among its energy density perturbation and 
the magnetic energy density: therefore, the source of (^]^) is now given by Ob which is 
no longer compensated, and the curvature grows. Comparing the above equation with 
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Eq. (86) of [^], we see that the mode sourced by the magnetic anisotropic stress is very well 
approximated by our analytical fit: the only difference among the two equations is the term 
—1/2 in Eq. (6.4) which becomes 1 — b/{^Ry) ~ 0.56 in Eq. (86) of |p. In our approach, 
however, we do not distinguish among the 'passive' and 'compensated' modes (see |Q): 
therefore, in Eq. ( |6.4| ) we get an additional contribution in the curvature proportional to 

Note that the curvature evolves with time at next-to- leading order in xi ^ 1, but 
this does not affect curvature conservation at large scales. First of all, the magnetic field 



contribution at order ©(l/xf) remains zero, as given in Eq. ( 4.23 ). Moreover, at next-to- 
leading order, the effect of the magnetic field on the curvature spectrum is constant in xi: 
|Cb(^)P oc x\, since the magnetic field anisotropic stress spectrum and the magnetic field 
energy density spectrum are flat for a causally generated field (up to the damping scale ko, 
see Eqs. (5.23) and ( |6.16| )). On the contrary, the inflationary contribution to the curvature 
spectrum scales as |Cinf(^)P x^^, and is therefore the dominant contribution at large 



scales. Hence, as already seen in section |4^, the curvature is still conserved at large scales, 
even in the presence of a magnetic field. 

We can now use the solution for Q at order 0{x^ to compute $ and ^ at the same order. 
Using Einstein's equation (|2.19| ), Eq. ( 4.22 ) can be rewritten as an evolution equation for 



2 



-3w'Hr]i 



C ■ 



(6.5) 



Before neutrino decoupling, tt,^ = and the source-term at order 0{x^) is given only by 
the term proportional to C, (Hb contributes at order 0{x^'^)). In order to determine the 
magnetic contribution to the Bardeen potential we therefore solve the above equation 
with C^iuB < y < Uu) 8LS source (given in Eq. ( |6.2| ) in the limit y <^ 1). The initial condition 
for ( |6.5| ) in this case is found by imposing the continuity of at i/b, ^^(j/b) = 0. Again 
in the limit He < y < Vu ^ ^, the solution reads 



<J>^(y) 



O(x0) 
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1 — log 

9 3^ 
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Vb) 2y 



3 ^ 95 ^ 

— JIb H IIb + log 

16 144 ^ 



^)nB 

yBj 



(6.6) 

In this expansion we have kept the term proportional to JIb y, since it corresponds to the 



solution given in the analytical estimate of [21|, Eq. (3.20), where however the magnetic 
anisotropic stress has been neglected^. 

Solving for $ after neutrino decoupling is more involved. Indeed, looking at ( |6.5| ) we 
see that the neutrino anisotropic stress at order 0{xi) contributes to the source term at 
order 0{xi). However, we cannot determine analytically tt^ at order 0{x1): this would 
require to redo the calculation presented in appendix^ at the following order, and therefore 



In the same way, ignoring the free-streaming neutrinos and taking the hmit y 3> 1, we would get 



C"(y): 
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where again the term proportional to Qb in $^ corresponds to the solution given in Eq. (3.20) of 
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to solve the whole system of coupled differential equations for the metric and the individual 
fluid components. This cannot be done analytically and is beyond the scope of this paper. 
We can however find an approximate solution by taking for tTu the following Ansatz 



'^Ay)\o{xi) = [di^B + d2nB j{y-yu) , (6.7) 

with di and d2 two arbitrary constants that we need to determine. This Ansatz has been 
chosen because it is consistent with the initial conditions of H, ^. It allows us to 
determine the first part of the source in ( |6.5D as a function of the constants di and d2 ■ The 
remaining part of the source is given by the curvature perturbation at y > yu, which we 
have determined previously and written in Eqs. ( |6.3| ) and ( |6.4D respectively in the limits 
y <^ I and y ^ 1. 



Knowing the source, we can solve Eq. (6.5) at order 0(x5) and determine ^ at this 



as 

mag) 



order as a function of di and d2. Note that ^ in this case corresponds to + ^ 
i.e. the magnetic part plus the neutrino magnetic one. The constant of integration is fixed 
by imposing that <I> is continuous at y^, i.e. <&mag(y!^) — ^^iVi^) derived above. Once 
$ is known, using Einstein's equations we can then derive ^ , D and V after neutrino 
decoupling as functions of di and d2. The constants di and d2 can then be determined 
solving the system of equations for the individual fiuid components: in particular, the 
energy and momentum conservation equations for neutrinos. The derivation is presented 
in appendix ^ We find then, in the limit <C 1 

[6 -4R^ + 8R,y log {yB / y u)] '^B - S Qb Ru , .2 2 on 

My)\oi..r = 2RM5 + 4R^) - ' ^'-'^ 

while for $ = <I)B _^ ^mag we find, for y,, < y <C 1 

^( Ru + [4(2 + R,) + 4(5 + 2R,) \og{y,/yB)] Hb 
= 2(15 + 4i?.) ' 



and for y 3> 1 

^iy)\o(xi) = -^^B + 



3 3 (yy 
— + - log — 

10 5 vyB 



Hb . (6.10) 



The full time dependence of the magnetic Bardeen potential ^(y) is rather complicated, so 
we do not write it here. It is shown in Fig. ^ where the transition from the radiation to the 
matter era is apparent. Note that, to solve both for the curvature and the Bardeen potential 
we made use of our fit for the neutrino anisotropic stress at leading order in krji <^ 1, see 
Fig. |l[ This allows us to find the complete solutions for the metric perturbations across 
the radiation-matter transition; however, it also causes minor differences among the early 
time limit of our analytical solutions and the corresponding initial conditions given in Q 
(for example, Eq. ( |6.9D slightly differs in the R^ dependence from the one that we could 
derive combining the passive and active modes given in the appendix of [^). 

From the above solution for $ we can compute ^, and in turns the Sachs Wolfe at order 
0{x^). This involves the computation of all the terms of Eq. (5.12), and its expression as 




Figure 2: Solid, blue: the magnetic contribution to the Bardeen potential at next-to- leading order 
^(y)\o{x^)' across the radiation-matter transition, together with its approximations at early and 



late time: Eq. (|6.9D, dotted, green and Eq. (3.10), dashed, red. In the plot, the magnetic energy 



density and magnetic anisotropic stress have comparable amplitude, and we have normalised the 
amplitude of <E> to IIb . 



a function of y is again rather complicated. We write it here for y = y^^c — 1 

1 

Bo(x;) ~ 540(15 + 4i?,)(fl,-l) 
+ 



AT 



495 + 2i?j,(811 + 224i?^ 



2960 - 2RJ977 + 108i?^) + 4(i?^ - 1)(505 + 108R^) log ( — 



(6.11) 



~ -0.2 JIb - 2.7Ub , 



where in the last line we have set Ru = 0.4, yi, = 10~^ and ye = 10^^^. Since the 
temperature anisotropy in the above equation is independent of k, it gives rise to a CMB 
spectrum scaling as £{£ + 1) Cf oc as we now demonstrate. 
The equivalent of Eq. (|5.20[ ) in this case is 



2i + l 



-0.2nB-2.7UB]je{kr]o 



(6.12) 



and the magnetic CMB spectrum becomes then (cf. Eq. j ^^M )) 



vr 



dkk'^ [0.04 |J^b(A:)P + 7.29 |nB(fc)P + 0.54 |CB(fc)|2] j|(A;r/o) , (6.13) 



where |nB(A:)p is given in Eq. ( [5.23 ), |ilB(^)P denotes the spectral amplitude of the mag- 
netic field energy density, 



(fiB(k)J7^(q)) = (2^)35(k-q)|J^B(A^)|' , 
and |Cb(A:)P denotes the spectral amplitude of the cross-correlation 

(f)B(k)n*(q)) = (2^)35(k-q)|CB(A:)P . 



(6.14) 
(6.15) 



in 

o 
o 

o 
I— ^ 
o 

o 
o 
o 
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The energy density spectrum and the cross-correlation one have been calculated in |l|, 
and they all share the same /c-dependence: similarly to Eq. (5.23), we set then 



in^m^ = n ^ ^ |Cb(A:)P = c ^ ^ 



Prad k 



(6.16) 



D 



Substituting the above equations in Eq. ( 6.13| ), we find 



£{i + l)Cf ~ -g^ivvcc) ro.04f] + 7.29 n + 0.54 C7l 
vr 



{B^f £{i + 1) rvoko 



Pri^d imkD f Jo 



dx X (x) 



g\Vrcc) [0.04 + 7.29 n + 0.54 C] T t 



(6.17) 



where for the last equality we have used approximation (A2) of Therefore, we confirm 
that the next-to-leading order contribution to the Sachs Wolfe effect behaves as found 



7. Conclusions 

In this work we present an analytical computation of the Sachs Wolfe effect induced by a 
primordial magnetic field. We have restricted our analysis to a magnetic field generated by 
a causal process, which acts on a time scale shorter than the Hubble time at generation, 
such as a first order phase transition in the early universe. The reason is that, under 
this hypothesis, the initial conditions for the metric and fluid variables are determined 
unambiguously by imposing conservation of the total energy momentum tensor across the 
transition, just as in the topological defects case []l3| , |l^ ]. After setting the initial conditions 
with this matching, we have solved Einstein's and conservation equations analytically and 
computed the magnetic effect on the gauge invariant metric and fluid perturbations in the 
large scales limit krji <^ 1, i.e. for scales larger than the horizon at recombination. Using 
these solutions, we have then determined the Sachs Wolfe effect. This is sourced only by the 
magnetic anisotropic stress: the temperature perturbation behaves as AT/T oc liB/ikrjiY 
at leading order in the krji <^ 1 expansion. A scaling of the CMB spectrum as £{£ + l)Cf oc 
£~^ trivially follows. 

However, as already pointed out in ||6|, |9|, the magnetic anisotropic stress is quickly 
compensated by the neutrino one, once neutrinos decouple from the primordial fluid. We 
have verified this compensation analytically, by solving the neutrino evolution equations in 
the limit kr]i <C 1, where we can consistently neglect multipoles higher than the anisotropic 
stress. Consequently, the contribution to the Sachs Wolfe at leading order in kr]i <C 1 gets 
strongly suppressed after neutrino decoupling, and the next-to-leading order in kr]i <C 1 
becomes the relevant one: the one which could in principle provide an observable effect 
in the CMB. In order to derive the magnetic Sachs Wolfe effect at next-to- leading order 
in kiji ^ 1 we had to compute the metric and fluid variables consistently at this or- 
der. This can be done analytically by solving the evolution equation for the curvature 
perturbation, and leads to a temperature anisotropy which is sourced both by the mag- 
netic anisotropic stress and the magnetic energy density: we found the approximate result 



AT/T ~ — 0.2r2B — 2.7Ii-Q. The new scahng of the CMB spectrum, arising from this 
contribution, takes the form (.{I + 1)C^ oc l"^. 

The matching conditions which we used to determine the initial conditions for the 
magnetic field perturbations imply that the curvature perturbation is continuous across 
magnetic field generation time. The subsequent evolution of the curvature perturbation is 
then completely determined. We found that, at leading order in krji <^ 1, the magnetic field 
does not contribute to the curvature, which is therefore conserved at large scales and given 
only by the inflationary mode (for which we have imposed adiabatic initial conditions). 
At next-to- leading order in kr]i <^ 1, the curvature perturbation is instead affected by 
the presence of the magnetic field. As already found in we confirm that, between 
the magnetic field generation time and neutrino decoupling, the curvature is sourced by 
the magnetic anisotropic stress and therefore grows (starting from zero at magnetic field 
generation time). After neutrino decoupling, when the magnetic anisotropic stress gets 
compensated, the part sourced by it remains constant, and the curvature perturbation 
gets a contribution from the magnetic energy density which is no longer compensated. 

Our analysis provides a clarification of the results present in the literature on the 
Sachs Wolfe effect arising from a primordial magnetic field. The magnetically induced 
CMB anisotropics have been studied in details in the literature by numerical integration 
of Boltzmann codes 0, ^, 0], and in most of the cases the large scale CMB spectrum 
behaves as £{i + l)Cf oc corresponding to a temperature anisotropy which does not 
depend on wavenumber. This is a consequence of the choice of the initial conditions. 
A Boltzmann code needs the initial conditions to be specified after neutrino decoupling, 
so the contribution to the Sachs Wolfe effect arising from the period of time between the 
magnetic field generation and the neutrino decoupling has been neglected in ||, ^, and 
the magnetic field anisotropic stress is exactly compensated by the free-streaming neutrino 
one since the beginning. This sets to zero the temperature anisotropy coming from the 
metric and fiuid perturbations at leading order in krji ^ 1 to start with. This choice of the 
initial conditions has been motivated by an analogy with the infiationary case, for which 
one can identify a growing and a decaying mode, and consistently neglect the decaying one. 
In fact, the infiationary solution for the Bardeen potential is, for y <^ 1 

J_ 6^ 20 oi ^ 3 r?i 20 
~ Ay^ xl 3 xl 4 r/ (krj)^ ^ 3 \r]J [ki])^ ' ^ ' 

where the mode proportional to bi is decaying with time, while the one proportional to 
ai is constant. One expects that from inflation the two modes are generated with the 
same amplitude: consequently, bi must be much smaller than ai since, at very early times 
{e.g. reheating), r]i/r]^ch ^ 1 [l^- The Sachs Wolfe effect from both modes at leading order 
in kr]i ^ 1 is given by (c/. also Eq. ( ^.10|) ) 

AT 



T 



12(1 + y^ec) hi 2(40 + 40y,ec + l^cc + 3j/D «i .72^ 
'yU2 + y,,,Yxl 3(2 + y,--^3 



The two modes have the same dependence on fe, and about the same amplitude for yrec — 1: 
therefore, knowing that 61 <C ai, one can neglect the contribution of the decaying mode 
(for an analysis of the effect of the infiationary decaying mode on the CMB, see p^ ). 



The situation for the metric perturbation induced by the magnetic field is a bit differ- 
ent. In this case, the solution at leading and next-to-leading order in krji <C 1 reads, for 
y <C1 (c/. also Eq. ( [Ol ) and Eq. Q) 



yZ 

= -3 



1 _ ^ 

y 

_ m 

V 



Hb 



+ 



+ 



'9 + 3^°K^J+^ 



5 2, 

1 — log 

9 3 ^ 



+ 



2ri 



Hb 



Hb 



(7.3) 



The term at leading order in krji <C 1 is decaying with time as y^^, whereas the next-to- 
leading order is growing logarithmically. One could therefore be tempted to neglect the 
first term on the right hand-side of the above equation. However, the key point here is that, 
contrary to the inflationary case, the two modes not only have no independent amplitudes 
(they are both proportional to Hb), but also do not share the same /c-dependence. Conse- 
quently, we see that the 'decaying' mode always has the same amplitude of the 'growing' 
one at horizon crossing A:t/ ~ 1, proportional to a few Hb- Since the CMB in general 
selects the contribution of modes which cross the horizon at recombination time r/i (the 
first peak), this means that the two modes can have the same impact on the CMB for 
kr]i ~ 1. Moreover, for the Sachs Wolfe effect, which selects modes still outside the horizon 
at recombination, we see that the 'decaying' mode could even dominate. 

However, the analyses of |0, ^, ^, do get the correct result: as we have seen, this 
is due to the compensating effect of the neutrinos, which exactly cancels Hb at leading 
order in krji <^ 1, strongly suppressing the amplitude of the 'decaying' mode. The relevant 
contribution to the CMB anisotropics is therefore caused by the 'growing' mode^. The 
correct CMB spectrum at large scales from the Sachs Wolfe effect from a causally generated 
magnetic field with spectral index n > 2 is therefore i{£ + ^)C^f oc Note however, that 
in order for this effect to be observable, the magnetic field should have an amplitude of 
at least a few nanoGauss: causal magnetic fields generated before Nucleosynthesis are 
unfortunately strongly constrained to amplitudes far too small to leave any observable 



effect in the CMB [23| 
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A. The metric and fluid perturbations at leading order 

In this appendix we present the solutions for the metric perturbations $ and ^, and the 
velocity V before and after the magnetic field generation, which are calculated via the 
matching procedure explained in section |^. 



^Note that, at next-to-leading order in k-qi <C 1 the cancellation from the neutrinos is not active: cf. the 
different amplitudes of Eq. (4.16) and Eq. (6.8) 



Before neutrino decoupling, from the expression of the density variable given in 
Eq. (4.11) one gets the Bardeen potential $ from Einstein's equation ( |2.18| ) 



= + *Bfe) = 2(8 + 3.)(10 +10. + V) a _ 24(1+ . to - .B) He 

The Bardeen potential ^' is obtained from Einstein's equation ( p.l9| ) 

= + * = = ^m±im±}^i±m 51 + 12te^-%.B -2.B) IlB , 

3{2 + yy xf y^{2 + yy xf 

(A.2) 

Note that the inflationary part is just opposite to The velocity perturbation V is 

obtained from Einstein's equation ( |2.20| ) 

V*iy) = V-jy) + V^(y) = + + '^'^ - ' HiL . (A,3) 

3(2 + !/)' II !/2(2 + !/)2 II 

After neutrino decoupling, to the above expressions one needs to add the neutrino 
contribution generated by the neutrino anisotropic stress. This can be calculated again 
from Einstein's equations starting from the additional contribution to the density variable 
given in Eq. ( [1.17 ). One gets then 



^ imR, {1 + y){2y + y,){y - y,f ai 24{1 + y){y - y,f Ub 

3(15 + 4i?^) y3(2 + y)3 3-2 y4(2 + y)3 xj ' 

^as ^ ^as^ ^ ^as ^^ 5^ 

^ 160 i^, yH4 + y)-3y^y,-il + y)yt ai 12{y - y,){y^ - 4y, - 3yy^) Ub 

3(15 + 4i?^) y3(2 + y)3 ^2 y4(2 + y)3 ' 

^as ^ -^as ^ p.|s 

1607?, (2y + y,)(y-y,)2 ai ^ 12(y-y,)2 Hb 



3(15 + 4i?,) y2(2 + y)2 y3(2 + y)2 ^.^ 

The above expressions depend on the particular form of the time dependence of the neutrino 
anisotropic stress that we choose, given in Eq. ( p^ . 



B. Resolution for the neutrino anisotropic stress at leading order 

We present here the derivation of the neutrino anisotropic stress at leading order 0{xi). 
We follow the method of In order to proceed analytically, we solve in the radiation 
era. The anisotropic stress so derived will be used to calculate the Sachs Wolfe effect 
and consequently it is assumed to be valid up to recombination time. Following Q, we 
derive a fourth-order differential equation for tTi, by combining Einstein's equations and the 
Boltzmann hierarchy 

D',,{y) + ^V,{y) = -W{y), (B.l) 

-Vliy) + l^y) - Dsuiy) = 4^(2/) , (B.2) 

xi 3 



e-i 



for 



>2 , 



(B.3) 
(B.4) 



2 

jH^in[ot(y) + + 3^^(u; - c^) - ^1 uMy)] , (b.5) 



where Htot = Ru'^u + SHb- We differentiate Eq. (|B.3|) and combine it with Eqs. ( [B.lD to 
(B.4) and with Einstein's equation ( 2.19D , in order to express and <I>" as functions of 



TTi, and its derivatives up to fourth order (they would be functions also of Dg y and higher 
multipoles, but these terms can be neglected since they are at higher order in xi <SC 1). 
<I>, and ^" so derived are then inserted back into Bardeen's equation ( [B.5D : at leading 
order, this equation becomes then a homogeneous fourth-order differential equation for the 
variable Htot 

Sy^ng; + 20y3ngl + 8i?,y2n;;^ _ 16i?,yCt + 16ii,ntot = . (B.6) 
Solving the above equation, we get 

/ N 3 ^ cos(Alogy) smiAXogy) n _ „, 

Ru ^/y y/y 



where A = \\J'^Ru — 1- Apart from the terms proportional to C3 and C4, this solution 
reproduces the one obtained in Eq. (16) of |^], which has been derived from a second order 
differential equation for t:^: Eq. (13) of the same reference. The time dependence of the 
terms proportional to C3 and C4 would not satisfy the second order differential equation 
given in P], which therefore does not admit, strangely enough, the standard inflationary 



solution for VTjy given for example in |16|. 



We determine then the constants ci to C4 by imposing the continuity of T^u^Vy^Dgy 
and C, at decoupling time yy, i.e. 



vr.(y.) = , (B.8) 
"5" 



<{y.) = -^Vy{yy) , (B.9) 



C(y.) = ^ , (B.ii) 



where for '^{yu) we take our solution '^^{y^) given in Eq. (|A.2| ) and evaluated well in 
the radiation era for y^ <C 1. Moreover, since we consider adiabatic initial conditions, at 
leading order in krji ^ 1 we have = V and Dgu = D. We take therefore for Vu{yv) and 



DsuiUu) our solutions V'^{yp) and D~^(y,y) given in Eqs. and (U). We then find for 

the neutrino anisotropic stress 



R,, 



40aiy^ 



+ 



+ 



15 + ARy 

3[16i;,(yB/y.-l) + 5]nB 
IQAR,, 



SHb , 40aiy^ 



2 n 



+ 



+ 



i?^ 15 + 4i?^ 
4ai(15 - 16i?; 



3(15 + 4i?^)y4 



//II y 

— COS A log — 

y \ Vu 



Vu ■ ( A, y 

— sm I A log — 



. (B.12) 



This solution is plotted in Fig. [l|, section 4^. We see that the neutrino anisotropic stress 
quickly adjusts to the opposite of the magnetic field one, and cancels it. It then starts to 
grow again, following the standard behaviour ~ 77^ coming from the inflationary solution 
and given for example in ||l^ . 

In order to solve the second order equation ( [4.1| ) for D analytically after neutrino 
decoupling, we introduce an approximated function which fits the behaviour in time of 
Eq. (IB.I2D , given by 



f ^ SHb fyl 
This is also plotted in Fig. 1^ section 4.2. 



40aiy(y - yu) 
15 + 4i?^ 



(B.13) 



C. Free-streaming neutrino contribution to the Sachs Wolfe plateau 



We evaluate here the contribution of free-streaming neutrinos to the Sachs Wolfe inflation- 
ary plateau, due only to their non-zero anisotropic stress (setting the magnetic field to 
zero). In the standard inflationary case this contribution adds to the purely inflationary 
result given in Eq. (|5.10| ). In Eq. (^.g|), it is represented by the part of — which is 
proportional to oi, as given in Eqs. ( [A.4[ ) and ( |A.5| ). In order to evaluate the Sachs Wolfe, 

) , as shown in Eq. (^.41) 



we simply have to differentiate Eq. (5 
expansion ©(l/xf), we find 



At leading order in the xi 



AT 



^inf 



^inf 



inf u 



160(yrec - Vu) 



2/rec 



2y3 - 2yrec(l + yu){yTcc + yu) 



Ru 0-1 



3yrec(2 + yr 



15 + AR^ x{ 



(C.l) 



Setting yy ~ 10 ^, yrcc — 1, -Ri/ — 0.4 we find 



AT 
AT 



inf 



u inf 



2.5 ai 



0.05 ai 



(C.2) 
(C.3) 



so that the neutrinos reduce the inflationary Sachs Wolfe plateau by about two percent. 



D. The metric and fluid perturbations at next-to-leading order. 



In this appendix we present the derivation of the metric and fluid variables at next-to- 
leading order, which we compute by solving Eq. ( |6.5D . We split each variable into its 
leading order component, which we calculated in Sec. [4.2| and which we label here with a o 
subscript, and its next-to-leading order component, which we label with a i subscript. We 
use the following Ansatz for the neutrino anisotropic stress 



T^i^liy) = ( dl^B + d2llB ) {y - Vvfxl , 



(D.l) 



and in order to compute the metric and fluid variables we need to determine the constants 
di and d2- To proceed analytically, we restrict to the radiation dominated era y <C 1. We 
neglect the matter components, and we split the fluid variables into a photon component 
and a neutrino component 



D = {1- R^)D^ + R^D^ , 
V = {1- R,)V^ + R,V, . 

The system of equations for these variables at next-to-leading order is 



1 — Ru)D^i + RyDjji 



$1 iv) 

D'yiiy) - TirjiDyi 



-3W RvTTul , 

Xi Xi 

4 o 



4x 



7I ) 



4x1 



T/ 2^ 

8X1 



(D.2 
(D.3 



(D.4 

(D.5 

(D.6 

(D.7 
(D.8 
(D.9 
(D.IO 
(D.ll 
(D.12 
(D.13 



Solving Eq. ( |6.5D with the ansatz ( p. ID , we can express <I>i as a function of di and d2- 
Eq. ( D.13| ) and Eq. (|D.5D give V^i, respectively ^'i as functions of di and ^2- Inserting 
these expressions in Eq. (p. 8]) allows to find the constant ^2 



3-2R^ , 41og(|^) 3 + 2di(15 + 4i;^) Ob 

2(15 + 4i?^) 



+ 



i?^(15 + 4:Ru) 15 + 4i?, 



(D.14) 
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With this, we find then 



T^uiiy) 



2RM5 + 4:R^) 



2(3-2i?^)nB + 8i?^log( — ) 



Hb - 3i?^f^B 



(D.15) 



We see that d\ and ^2 are not independent and that ^2 determines completely the solution. 

We can then write the full solutions for the metric and fluid variables in the limit 
Ub < yu < y ^ ^, which would correspond, in our approach, to the initial conditions for 
the Boltzmann hierarchy given in [Q, ^, ^ in the syncronous gauge 
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